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Abstract. The Ising model is widely regarded as the most studied 
model of spin-systems in statistical physics. The focus of this paper is 
its dynamic (stochastic) version, the Glauber dynamics, introduced in 
1963 and by now the most popular means of sampling the Ising measure. 
Intensive study throughout the last three decades has yielded a rigorous 
understanding of the spectral-gap of the dynamics on Z 2 everywhere 
except at criticality. While the critical behavior of the Ising model has 
long been the focus for physicists, mathematicians have only recently 
developed an understanding of its critical geometry with the advent of 
SLE, CLE and new tools to study conformally invariant systems. 

A rich interplay exists between the static and dynamic models. At the 
static phase-transition for Ising, the dynamics is conjectured to undergo 
a critical slowdown: At high temperature the inverse-gap is O(l), at the 
critical /3 C it is polynomial in the side-length and at low temperature it 
is exponential in it. A seminal series of papers verified this on Z 2 except 
at /? = /3 C where the behavior remained a challenging open problem. 

Here we establish the first rigorous polynomial upper bound for the 
critical mixing, thus confirming the critical slowdown for the Ising model 
in Z 2 . Namely, we show that on a finite box with arbitrary (e.g. fixed, 
free, periodic) boundary conditions, the inverse-gap at ft = /3 C is poly- 
nomial in the side-length. The proof harnesses recent understanding of 
the scaling limit of critical Fortuin-Kasteleyn representation of the Ising 
model together with classical tools from the analysis of Markov chains. 



1. Introduction 

The classical Ising model on the lattice is one of the most studied models 
in mathematical physics with thousands of research papers since its intro- 



duction in 1925. In his famous work from 1944, Onsager 44 exactly solved 
the model in two dimensions, thereby determining its critical temperature. 
Ever since, physicists directed most of their attention to the fascinating be- 
havior of the model at criticality (see for instance the 20 volumes of |12|). 
In this regime the model on Z 2 exhibits delicate fractal geometry whose 
rigorous understanding was only recently obtained using conformally invari- 
ant scaling limits. Here we report the first rigorous confirmation that the 
spectral-gap of the Glauber dynamics for the critical Ising model in Z 2 , 
perhaps the most practiced methods for sampling its Gibbs distribution, is 
polynomial in the side-length. 
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The Glauber dynamics (also known as the stochastic Ising model) is a 



family of Markov chains introduced by Glauber 17 in 1963, which both 
models the dynamic evolution of the Ising model and provides a simple 
algorithm for sampling from its stationary distribution. The most notable 
examples are the heat-bath dynamics and metropolis dynamics, both highly 
used in practice thanks to their simple and natural transition-rules. 

An extensive program of work by mathematicians, physicists and com- 
puter scientists has related the static spatial-mixing properties of the Ising 
model to the mixing-rate of the Glauber dynamics, measured in terms of 
the gap in the spectrum of its generator. Relying on many experiments 
and studies in the theory of dynamical critical phenomena, the spectral-gap 
of the dynamics on a finite box in the lattice is conjectured to have the 



following critical slowdown behavior (e.g., [20 30 38 60|): 

• At high temperatures the inverse-spectral- gap is O(l) 

• At the critical /3 C it is polynomial in the side-length 

• At low temperatures it is exponential in the side-length 

(in dimensions d > 3 the surface-area plays the role of the side-length) . As 
we detail later, in a long series of seminal papers over the last three decades 
this rich interplay between the static and dynamic models was confirmed 
for the 2-dimensional Ising model with the crucial exception of the critical 
temperature (5 C . By contrast to the detailed picture by now known for 
(3 ^ (3 C , the behavior of the inverse-gap at criticality remained a stubborn 
and fundamental open problem, with no known sub-exponential bounds. 

Our main result establishes that the inverse-gap is indeed polynomial in 
the side-length, thereby confirming the critical slowdown behavior of the 
Ising model on W? . 

Theorem 1. Consider the critical Ising model on a finite box A C Z 2 of 
side-length n, i.e. at inverse-temperature j3 c =\ log(l+\/2)- Let gap^ be the 
spectral- gap in the generator of the corresponding continuous-time Glauber 
dynamics under an arbitrary fixed boundary condition r. There exists an 
absolute constant C > (independent of A, r) such that (gap^) -1 < n . 

Using well-known relations between L 1 and L 2 mixing (see, e.g., [48]) an 
immediate corollary of our main result is that the total-variation mixing 
time of the Glauber dynamics is bounded above by a polynomial in the 



side-length n (see Corollary 2.1 in Section 2.2) 



Furthermore, we show a generalized version of the above theorem (see 



Theorem 3.1) that applies to rectangles of unbounded aspect-ratio (e.g., an 
nxe" box), where the polynomial upper bound depends only on the shorter 
side-length. 
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(a) Critical Ising model 



(b) Critical FK-Ising 



Figure 1. The FK-representation for the critical Ising model. 

(A) Critical Ising on a 256 x 256 lattice; 64 x 64 box highlighted. 

(B) Coupled FK-configuration corresponding to highlighted box. 



Lower bounds on the inverse-gap can be derived already from the work 

demonstrating that it grows at least 



of Onsager, as observed by Holley 21 



polynomially fast with n. Following this approach we complement the above 
upper bound by showing that the inverse-gap is of order at least n 7//4 (see 



Theorem 4.1 ). 



Remark. No attempt was made to optimize the exponent C > given in 
Theorem[T} As we later elaborate, our proof determines this bound explicitly 
in terms of a bound c + < 1 on the probability of a certain crossing event 
in the critical Fortuin-Kasteleyn representation of the Ising model (namely, 
the dependence on 1/(1 — c + ) is logarithmic). The true exponent is believed 



to be universal and numerical experiments supporting this [ 18, 26, 43 ,60, 61 
(among others) suggest that its value is about 2.17. 



1.1. Background and previous work. In what follows we give a brief 
account of related work on the stochastic Ising model. For a more extensive 
exposition, see e.g. 34,38 . 

The classical Ising model on a finite box A in the square lattice Z 2 with 
no external field is defined as follows. Its set of possible configurations is 
Qjs = where each configuration corresponds to an assignment of 
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plus/minus spins to the sites in A. The probability that the system is in a 
configuration a 6 VL IS is given by the Gibbs distribution 



Ma(o-) = 7? exp I f3 V a(u)a{v) I , (1.1) 




where the partition function Z is a normalizing constant. The parameter /3 
is the inverse-temperature of the system; for (3 > we say that the model 
is ferromagnetic, otherwise it is anti-ferromagnetic. 

These definitions extend to the infinite volume lattice Z 2 , where there is 
a critical point ji c =\ log(l + \/2) such that the Gibbs distribution is unique 
if and only if /? < f3 c . We note that the focus of this work is restricted to the 
case of no external field (h = 0) since otherwise there is no phase-transition 
in the above setting and the dynamics for the system is well understood. 

Define the boundary of a set Ac V, denoted by dA, as the neighboring 
sites of A in V \ A and call r £ {±1}^ A a boundary condition. The Gibbs 
distribution conditioned on r will be denoted by and put A = AU<9A. A 
periodic boundary condition on a box corresponds to a 2-dimensional torus. 

The Glauber dynamics for the Ising model is a family of continuous-time 
Markov chains on the state space Q IS , reversible with respect to the Gibbs 
distribution, given by the generator 

(£f)(a) = Y J c(x,a)(f(a*)-f(v)) (1-2) 

xeA 

where a x is the configuration a with the spin at x flipped. The transition 
rates c(x, a) are chosen to satisfy finite range interactions, detailed balance, 
positivity and boundedness and translation invariance (see Section [2]). As 
mentioned above, the two most notable examples for the choice of transition 
rates are 

(i) Metropolis: c(x, a) = exp (^2f3a(x) Y^y~ x A 1. 

(ii) Heat-bath: c(x, a) = 1 + exp f — 2/3a(x) Yly~ x a (y) 

These chains have intuitive and useful graphical interpretations: for in- 
stance, the heat-bath Glauber dynamics is equivalent to updating the spins 
via i.i.d. rate-one Poisson clocks, each time resetting a spin and replacing it 
by a sample according to the conditional distribution given its neighbors. 

Perhaps the most fundamental property of the dynamics is the gap in 
the spectrum of its generator, which governs the rate of convergence of the 
Glauber dynamics to equilibrium in L 2 (fi). The spectral-gap is defined via 
the following Dirichlet form: 

gaP = fvar(7)' (L3) 
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where the infimum is over all nonconstant / £ L 2 ([i) and 

8(f) = (£/, f) LHfl) = \Y, ti°)<*, *) [/("*) " ■ (1-4) 

Alternatively, one may define the gap via the discrete-time analogue of the 
chain, in which case it is given in terms of the largest nontrivial eigenvalue 
of the transition kernel. 

We next present a partial list of results obtained on the spectral-gap in 
various temperature regimes over the last half-century. 

High temperature regime. A series of breakthrough papers by Aizenman, 
Dobrushin, Holley, Shlosman, Stroock et al. (cf., e.g., (!} [TTJ [2lJ [23j [24| [34} 



35, 40 42 56 58 64, 65)) beginning from the late 1970's has developed the 
theory of the convergence rate of the Glauber dynamics to equilibrium. 

Aizenman and Holley [l] showed that the spectral gap of the dynamics on 
the infinite-volume lattice is uniformly bounded whenever the Dobrushin- 



Shlosman uniqueness condition holds. Stroock and Zegarlihski 56 , 58 64 
proved that the logartihmic-Sobolev constant is uniformly bounded given 
the Dobrushin-Shlosman mixing conditions (complete analyticity) . Finally, 



in 1994 Martinelli and Olivieri 40 , 41 extended this result to cubes under 



the more general condition of regular complete analyticity, shown to hold 
for any /3 < (3 C in two dimensions. In particular, this confirmed the con- 
jectured behavior of the inverse-gap in the high temperature regime in Z 2 . 
This spectral-gap was shown to determine the precise asmyptotics of the 



L 1 mixing time of the dynamics in the recent work 36 , where the authors 



established the cutoff phenomenon for the dynamics in this regime. 



See the excellent surveys 38,39 for further details. 



Low temperature regime. In this regime, in accordance with the conjectured 
critical slowdown behavior of the Ising model, the Glauber dynamics under 
the free boundary condition was expected to converge exponentially slowly 
in the side-length. This was first established by Schonmann [49] in 1987 
for sufficiently low temperatures as a direct corollary from a large deviation 
estimate on the magnetization in a square. This large deviation result was 
extended to any (3 > f3 c by Chayes, Chayes and Schonmann [7j, implying 
that the inverse-gap under free boundary conditions is at least exp(— cn) for 
some c = c(/3) > 0. This result, concurring with the projected behavior of 
the inverse-gap in this regime, also appeared explicitly and independently 
by Thomas [59] . 

In 1995, precise large deviation rate functions for the magnetization were 
established by Ioffe |25| for all /3 > /3 C . Building on the work of Martinelli |37| 
(treating low enough temperatures), this then culminated in remarkably 
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sharp estimates for the inverse-gap throughout the low temperature regime 
in the work of Cesi, Guadagni, Martinelli and Schonmann |6j in 1996. 

Critical temperature. Evidence that the inverse-gap is at least polynomial 
in the side-length at the critical temperature followed from the polynomial 
decay in the spin-spin correlation, whose asymptotics were determined by 
Onsager [44 1. On the other hand, despite the vast body of papers that 
accumulated on the non-critical regimes, no sub-exponential upper bound 
was known for the inverse-gap at criticality. 

In fact, the only geometries where the inverse- gap of dynamics for the 
critical Ising model was shown to be polynomial (in the appropriate param- 
eter for the given underlying geometry) were the complete graph j9], where 
the large symmetry renders the methods of proof useless for lattices, and 



the regular tree 10 , whose non- amenability and non-transitivity also re- 
quire geometry-specific methods. Even for these much simpler geometries 
the proofs were highly nontrivial. 

Indeed, the complexity of the critical behavior can be witnessed by the 
scaling limit of the static Ising model. The understanding of the limit 
emerged in the last decade, pioneered by the introduction of the Schramm- 
Loewner Evolution (SLE) by Schramm [50] (see |62| for more on SLEs). This 
powerful machinery has revolutionized the study of critical phenomena in 
two dimensions and allowed direct calculation of various critical exponents 



(see for instance 31 , 32 , 54 to name just a few). 

Striking new results by Smirnov (cf. |53,55|) describe the full scaling limit 
of the Ising cluster interfaces at criticality as the Conformal Loop Ensemble 



(CLE) with parameter k = 3 (also see [33 52 63] ). 

An important role in the development of the scaling limit theory for the 
Ising model was played by its counterpart, the Fortuin-Kasteleyn representa- 
tion (formally defined in Section [2]). Scaling limit results that were initially 
obtained for this model at criticality were thereafter translated to the corre- 
sponding Ising measure. For instance, the full ensemble of cluster interfaces 
in the FK-Ising model converges to a (nested) CLE K for k = 16/3. 

Along side the ground breaking results on critical percolation models in 
two dimensions, new tools were developed to study conformally invariant 
systems. A recent application of this theory yielded Russo- Seymour- Welsh 
type estimates ( [47|[5l] ) for the crossing probability in a rectangle of bounded 
aspect ratio in the critical FK-Ising model under arbitrary boundary con- 
ditions, recently obtained by Duminil-Copin, Hongler and Nolin [13) . See 
also [8] for this result under a specific boundary condition, as well as [5] 
for an argument inferring the general estimate from announced results on 
the convergence of the scaling limit of the spin cluster boundaries to CLE3 
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and from its Brownian loop soup representation. As we later state, the 
RSW-estimate of [13] is a key ingredient in our proof. 

As a consequence of Theorem[TJ which verifies the conjectured polynomial 
behavior of the inverse-gap at criticality, the critical slowdown of the Ising 
model in 2 dimensions is now fully established. 

1.2. Polynomial mixing time under arbitrary boundary conditions. 

Recall that Theorem [T] bounded the inverse-gap for any fixed boundary 
condition. Our proof in fact extends to free or periodic boundary conditions 



as well as mixed boundary conditions (see Theorem 4.2). Furthermore, by 



a standard reduction the result also carries to the anti-ferromagnetic Ising 



model (see Corollary 4.5) 



To the best of our knowledge, this gives the first rigorous polynomial- 
time algorithm for approximately sampling the critical Ising model in Z 2 , 
also leading to an approximation of its partition-function, under arbitrary 
(e.g. mixed) boundary conditions. (In the absence of boundary conditions 
Jerrum and Sinclair provided an efficient approximation scheme for the Ising 



partition-function on any graph in their celebrated work 27 . Based on this 
algorithm Randall and Wilson 1461 gave an efficient approximate sampler for 
the Ising model on general graphs in 1999, applicable whenever the boundary 
conditions are free/all- plus/all- minus or periodic.) 

1.3. Efficient rigorous perfect simulation. In their highly influential 
work |45j, Propp and Wilson introduced the method Perfect Simulation (or 
Coupling From The Past) to exactly sample from the stationary measure of 
certain Markov chains. Perhaps their most prominent application for this 
method was sampling the Ising model at criticality via Glauber dynamics. 

While the Propp- Wilson algorithm is guaranteed to produce a precise 
sample from the Ising measure, it was not rigorously shown to be efficient, 
that is to have an expected running time polynomial in the size of the 
lattice. The only guarantee given for the running time was in terms of 
the total-variation mixing time, which was unknown at criticality. Note 
that experimental results found the method to be extremely effective in 
this regime (a fact that served as additional supporting evidence that the 
inverse-gap of the dynamics is polynomial at j3 = j3 c ) . 



By the new results in this work (namely, Corollary 2.1) we now have 
the first proof that the Propp- Wilson algorithm runs in polynomial time. 
This constitutes the first rigorously proven efficient algorithm (in addition 
to being very simple) for perfectly simulating the critical Ising model on 
the square lattice. Furthermore, this algorithm is valid under arbitrary 
boundary conditions. 
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1.4. Main techniques. A common feature of most analyses of the Glauber 
dynamics is to utilize the spatial-mixing properties of the static Ising model 
in order to control the mixing rate of the dynamical process. In the high 
temperature regime, this is typically done by measuring the influence of indi- 
vidual boundary condition spins on sites and its rate of decay with distance. 
However, at the critical temperature, the slow decay of correlations precludes 
this approach. Indeed, as demonstrated by the spin-spin correlation result 
of Onsager and further illustrated by the existence with positive probability 
of "large" conformal loops as Ising cluster interfaces, there are long-range 
correlations at criticality which foil the standard coupling techniques. 

Instead, we apply ideas from the study of the conformal invariance of 
the Ising model, and crucially the RSW-type estimate of 



13 



to obtain the 

spatial-mixing result required for our analysis (Theorem 3.2). Rather than 



considering the effect of a single spin on the boundary, here we analyze the 
effect of an entire face of the boundary, deducing just enough spatial mixing 
to push our program through, with the help of additional ingredients from 



the analysis of Markov chains, in particular the block-dynamics method 38 



2. Preliminaries 

Throughout the paper we will use the notation [x, y\ = [x, y] nZ, whereby 
for instance a square lattice of side-length n can be denoted by [[l,n] 2 . 

2.1. Glauber dynamics for the Ising model. The Glauber dynamic for 
the Ising model on a finite box A C Z 2 , whose generator is given in (1.2), 



accepts any choice of transition rates c(x, a) which satisfy the following: 

(1) Finite range interactions: For some fixed R > and any x G A, if cr, u' G 
Ois agree on the ball of diameter R about x then c(x, a) = c(x, a'). 

(2) Detailed balance: For all a G fiig and x G A, 

c(x, a) 



(3) Positivity and boundedness: The rates c(x, a) are uniformly bounded 
from below and above by some fixed C\, C2 > 0. 

(4) Translation invariance: If a = o~'(- + £), where £ £ A and addition is 
according to the lattice metric, then c(x, a) = c(x + I, a') for all x G A. 

The Glauber dynamics generator with such rates defines a unique Markov 
process, reversible with respect to the Gibbs measure n\. 
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2.2. Mixing in L 1 and in L 2 . For any two distributions 99,^ on Q, the 

total-variation distance of (p and ifi is defined as 



sup \ip(A) - $(A)\ = - J2 - i>(x)\ ■ 



\<P - tI>\\tV = - . ,r ■ 



The (worst-case) total-variation mixing-time of an ergodic Markov chain 
(X t ) with stationary distribution fi, denoted by t Mlx = t Mlx (l/e), is 

min < t : max ||Pa.(Xt 6 •) — /u||tv < V e f > 

where P x denotes the probability given that Xq = x. It is easy and well 
known (cf., e.g., 48 ) that t mx < gap -1 log where ^ min = mm xen K x )- 
Together with Theorem [TJ this immediately implies the following: 

Corollary 2.1. Consider the critical Ising model on a finite box AcZ 2 of 
side-length n, i.e. at inverse-temperature (3 C = 2 l°g(l + v2). Let i^ IX 6e i/ie 
worst-case total-variation mixing time of the corresponding continuous-time 
Glauber dynamics under an arbitrary fixed boundary condition r. There 
exists an absolute constant C > (independent ofA,r) such thatt^^ < n c . 

2.3. Single-site vs. Block dynamics. Let B = {B\, . . . , B^} where the 
Bi's are subsets of sites to be referred to as "blocks". The (continuous- 
time) block dynamics corresponding to B is the following Markov chain: 
Each block corresponds to a rate-one Poisson clock and upon it ringing we 
update the block according to the stationary distribution given the rest of 
the system. That is, the entire set of spins of the chosen block is updated 
simultaneously, whereas all other spins remain unchanged. One can verify 
that this dynamics is reversible with respect to the Gibbs distribution. 

The following proposition reveals the remarkable connection between the 
single-site dynamics and the block dynamics. 

Proposition 2.2 ( |38[ Proposition 3.4]). Consider the Glauber dynamics 
for the Ising model on A C Z 2 with boundary condition r. Let gap A be the 
spectral-gap of the single-site dynamics on A and gap^ be the spectral-gap 
of the block dynamics corresponding to B\, . . . , B^, an arbitrary cover of A. 
The following holds: 

where N x = #{i : Bi B x}. In particular, 

gapX > ( sup #{i : B { 3 x}) gap^ inf inf gap£. . 
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A notable example for key properties of the Glauber dynamics that were 
established using this approach is the estimate on the spectral-gap through- 
out the high temperature regime (cf., e.g., [38]). 

2.4. FK model. Invented by Fortuin and Kasteleyn [15] around 1969, the 
FK-model (also known as the random- cluster model) with parameter q for 
a graph A with edge set E is defined as follows. Its state space is all assign- 
ments of open/closed values to the edges (resp. 1 and values), where the 
probability of a configuration oj 6 {0, 1} E is given by the FK-measure: 

utSw) = i p °H(i -p) c H g fc H , 

where o(oj),c(uj) denote the number of open and closed edges in lo resp., 
k(uj) counts the number of clusters and Z is a normalizing constant. 

The case q = 2 is closely related to the ferromagnetic Ising measure, as 



detailed by the following coupling due to Edwards and Sokal 14 , hence 
this special case is also known as FK-Ising (or the FK-representation of the 
Ising model). Let FK and fi IS denote the state spaces of the Ising and 
FK measures respectively; the joint probability of a £ J2 IS and oj 6 PK 
according to the coupling is 

F(<J,0j) OC Y[ [(1 -P)l{w(e)=0} +P 1 Me)=l} 1 M:r)= ( T(y)}] , 

where the relation between < p < 1 and /3 > is given by 

p = 1 - exp(-2/3) 

(note that this sometimes appears as p = 1 — exp(— /3), e.g. in the context of 
Potts measures, resulting from a slightly different normalization of the Ising 
Hamiltonian). Unless stated otherwise, set (3 and p to their critical values 

P = Pc = - log(l + V2) , p = p c = Psd = 1 + ^ ' 

A useful corollary of the coupling is that an Ising configuration a ~ 
can be obtained from the FK-Ising representation oj ~ by assigning i.i.d. 
uniform spins to the clusters. 

One can also consider the FK-measure under boundary conditions. Here 
the constraints are "wiring" of sites, thereby affecting the cluster-structure of 
oj (wired sites are part of the same cluster) . Given a wiring r we denote by v\ 
the FK-measure with respect to this boundary condition. The two extreme 
boundary conditions are the "wired" r = 1 (all boundary vertices are pair- 
wise connected) and the "free" r = (no wiring) boundary conditions. 
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The FK- model is monotone with respect to boundary conditions: In the 
natural partial order on the wirings, whenever £ < rj the measure v\ stochas- 
tically dominates lA. This property is particularly useful when combined 
with the so-called Domain Markov property, asserting that ^ (w 6 • 
for some subset of the edges A C E is precisely the FK-measure on the re- 
maining (unexposed) edges under the boundary condition induced by the 
wirings of £ and a; (A). 

A final ingredient which we will need is the FKG-inequality, stated next 
for the special case of the FK-measure: 

Theorem 2.3 ([16,22 , special case). For any graph A and boundary con- 
dition £, if X,Y : f2 FK —> M are two increasing (decreasing) functions then 

4(XY)>u{(X)4(Y). 

Throughout the paper we will use the term cluster for its usual interpre- 
tation as a maximal subset of sites connected by open paths, including via 
the boundary wiring. As a distinction, we will use the term component to 
refer to a maximal subset of sites that are connected by open paths without 
using the boundary wiring. 

See e.g. 19 and [3] for further information on the FK-model. 



3. Polynomial mixing at criticality 

In this section we provide the proof of Theorem [TJ establishing that the 
inverse-gap of the Glauber dynamics for the critical Ising model on a finite 
square lattice is polynomial in the side-length. We will in fact prove the 
following stronger version of Theorem [T] which allows for rectangles with 
unbounded aspect-ratio and bounds the inverse-gap solely in terms of their 
shorter side-length. 

Theorem 3.1. Consider the critical Ising model on a box A C 1? with 
dimensions m x n. Let gap^ be the spectral- gap in the generator of the 
corresponding continuous-time Glauber dynamics under an arbitrary fixed 
boundary condition r. There exists an absolute constant C > (independent 
of A, t) such that for any m = m(n) we have (gap^)" 1 < n c . 

For instance, the inverse-gap for the critical Ising model on an extremely 
long rectangle, e.g. an n x exp(exp(n)) box, is bounded by the exact same 
n bound given for the square. 



3.1. A spatial-mixing result and proof of Theorem 3.1 We derive 



our bound on the inverse-gap from the following spatial-mixing result which 



is of independent interest. Its proof is given in Section 3.2 
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Theorem 3.2. Let A = [l,r] x [l,r'J for some integers r,r' satisfying 
r'/r > a > with a fixed and let A T = [l,r] x [pr, r'J for some p satisfying 
a < p < r'/r. Let £,77 be two boundary conditions on A that differ only on 
the bottom boundary [l,r] x {0}. T/ien 

Mi(«r(A T ) G ■) - Ma( ct ( A t) G •) < exp(-tfp) , 

TV 

where 5 > is a constant that depends only on a. 

Proof of Theorem [TJ We proceed by using the method of Block Dynamics 
to recursively relate the spectral-gaps of boxes of decreasing sizes, following 
the approach of Martinelli [38] . 

Let A = [l,r] x [l,r'] be a box of dimensions r x r' and without loss of 
generality assume r' > r. For some I G {1, . . . , \_-J r'/r\ } partition A into 
two vertically overlapping boxes: 

Ai=Ai(i) = [l,r] x [I r ' + ^±vV,r'] 
A 2 =A 2 (*) = [l,r] x [l^r' + fx/rV] 

(e.g., if r' = r then these overlap in the middle third of A). 

Denote by gapg the spectral-gap of the block-dynamics on A (see Section^ 
for the definition of this dynamics) corresponding to the blocks B = {Ai, A2} 
and the boundary condition £. 

Claim 3.3. For any boundary condition £ on the box A, the spectral-gap 
of the block- dynamics that corresponds to B = {Ai,A2} as defined above 
satisfies gap^ > 1 — exp(— cy/r'/r), where c > is an absolute constant. 

Proof. It is sufficent to consider the discrete-time version of the dynamics, 
whose spectral-gap we will denote by gap = ^gap^. Let (X t ) and (Y t ) be 
two instances of the discrete time block-dynamics on A under the boundary 
condition £, started at two arbitrary initial positions Xq and Yq respectively. 

Couple the choices of updated blocks in (Xt) and (Yt) and consider the 
case where two distinct blocks are updated at times t, £ + 1 for some t > 1. 
If we first update Ai then upon updating this block the boundary condition 
is the same given Xt—i, Yt—\ except (possibly) on the boundary with A \ Ai. 



Using a maximal coupling we have by Theorem 3.2 (setting parameters 



a = I and p = \sjr ! jr) that with probability at least 1 — exp(— bp) we 
can couple Xt and Yt so that they agree on [l,r] x [|r' + |y / rr', r'\. In 
particular they would agree on A\A 2 . On this event, upon updating A 2 in the 
following step the two chains may be coupled with probability 1. Similarly, 
if we first update A 2 then for the same a, p we have a probability of at least 
1— exp(— Sp) to couple X t , Y t so that they agree on [1, rj x [1, ^r' + e ^\/rr'\, 
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again implying coalescence of Xt+i, Yt+i. We deduce that for any t > 2 

nXt ?Y t )<(l-l(l- e~ s ") ) ^ = Q (l + e" W»>^ 
and in particular the exponential decay of the coupling time satisfies 

1, „™ ,„ , ™ w, 1 



lim ^logmax||P x (X t e -)-F y (X t e -)IItv < ^ (l+eW'/r . (3.1) 

As the exponential decay of (Xt) to equilibrium is governed by A, the second 
largest (in absolute value) eigenvalue of the transition kernel, we have that 
A is precisely the above limit, hence 

gap« = 2gap = 2(1 - A) > 1 - e ~^^~ r . 
This completes the claim. ■ 



To obtain Theorem 3.1 we now relate the spectral-gap of A with that of 
the Aj's. By the above claim and Proposition 2.2 we obtain that for some 
absolute c > and any choice of I we have 



(gap!)- < E ^ ( 4 f* (l - e-W^Tr)" 1 max ( gapLm )- 



1. 1 1 



where 4 = Z ff n{(p-)c(x,*)\f(o*) - f{a)f and N x {i) = #{* : M*) 3 4- 
Averaging over the L = [\/r' /r\ admissible values of £ now yields 

(gapi) < z 2^ — (1 - e V / j max(gap^ 



s( S E¥ i )( 1 -^)" I X^W- 1 



-1 



where the last inequality is due to the fact that (crucially) every x G A 
appears in the overlap Ai(£)nA2(^) for at most one value of £ Observe that 
since each Aj(£) has dimensions r x r" for some gr' < r" < |r' we can write 

(gapi)" 1 < (l + — (l-e-^^rnaxma^gapl,)- 1 , (3.2) 

where A' runs over all boxes of the form [1, r\ x [1, r"\ with < r" < |r'. 
Therefore, as long as r" > r we can recursively apply the above argument 
and obtain that 

(gapi) -1 < ^ maxmax(gap](,) _1 , 
A' V 
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where A' goes over all boxes of the form [l,r] x [l,r"]] with |r < r" < r 
and where 

oo 1 

^o = Il(l + (3/2)" fc / 2 ) (l - e - c ^ k/2 ) ' < oo . 

k=0 



An additional application of the recursion (3.2 ) (this time reversing the roles 
of r, r") now yields 

(gap!) -1 ^ i9maxmax(gap](,)" 1 , 
A' V 

where i? = 2i?o/(l — e_c ) an d A' runs over all boxes of the form [1, a] x [1, 6J 
with a A 6 < |r. 

The proof of Theorem |3.1| is now concluded by repeatedly applying this 
argument log 3 / 2 n times inductively, implying that the inverse-gap of the 
single site dynamics on an m x n square lattice with arbitrary boundary 
conditions is at most $ log3 / 2?1 = n log3 / 2,? , as required. ■ 



3.2. Proof of Theorem 3.2 (spatial mixing for critical Ising on Z ). 

The following is an immediate corollary of the Edwards-Sokal |14| coupling 
of the Ising and FK-Ising models under free boundary conditions: 

Fact 3.4. Let A C 1? be a finite box and let £ be a wiring of its boundary 
vertices dA. Let S\,S2,--- denote the clusters of dA that are induced by 
the boundary condition £. Let oj £ ri FK (A) be distributed according to the 
FK-Ising measure ui and transform it to a spin configuration a G IS (A) 
by selecting an i.i.d. spin per cluster. Then a is distributed according to ^ 
conditioned on every Si receiving the same spin in each of its sites. 

Indeed, to see the above simply observe that wiring two given sites u, v is 
equivalent via the Edwards-Sokal coupling to introducing an interaction of 
strength f3(u, v ) = oo between these sites in the corresponding Ising model, 
i.e. conditioning that they obtain the same spin (see Fig. [2]). 

An arbitrary boundary condition £ for the Ising model on A induces a 
partition of the vertices of dA into two subsets, V(£),M.(£), corresponding 



to the plus-spins and minus-spins respectively. With Fact 3.4 in mind, by 
slight abuse of notation we will also let £ denote the boundary condition for 
the FK-model generated by wiring all the vertices of V amongst themselves 
pair- wise and doing the same for those in M. Further define the following 
event for configurations U) G f2 FK (A): 

At = {r(OS>M(0} , 

that is, cj does not admit an open path that connects the two subsets V,M. 

4 



By the above fact, if we condition that a configuration oj ~ v\ satisfies the 
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(a) Ising model 



(b) FK-Ising model 



Figure 2. Coupling of the Ising model and its FK-Ising represen- 
tation under arbitrary fixed boundary. Blue sites and red sites are 
each wired pair- wise in the FK-configuration on the right. 

event A^ then the Edwards-Sokal procedure for generating an Ising config- 
uration from it agrees with the Ising boundary condition £ with probability 
2 independently of uj. We thus arrive at the following corollary: 



Corollary 3.5. Let £ be an arbitrary boundary condition for the Ising model 
on A and define and as above. Let u G PK be distributed 

according to i4 (• | A^) and produce an Ising configuration a 6 $7i S from uj 
by assigning a plus-spin to the cluster of V , a minus-spin to the cluster of 
Ai and i.i.d. uniform spins to all other clusters. Then a ~ [A. 

Note that for some boundary conditions £ the probability of A^ might be 
extremely small, e.g. exponentially small in the side-length of A, in light of 
which treating this conditional space is quite delicate. 



Let £, rj be the Ising boundary conditions given in Theorem 3.2 and define 

r^[0,r + l] x{0}, 
while recalling that by definition £, n differ at most on the sites of T. In light 



of Corollary 3.5, consider the FK-Ising model on A = [0, r + 1] x [0, r' + 



under the two wirings corresponding to £, rj as defined above (wiring the 
sites in V pairwise and the sites in M. pairwise). 

We next wish to construct a coupling of the measures 



{4(- | A*) , v\{- I A") , 4} , 



where is the FK-measure under the (fully) wired boundary condition. 
This will be achieved by gradually exposing all the sites of the configuration 
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(jl ~ u \ are connec ted by an open path to the bottom boundary T. 
Let E = E(u 1 ) denote this set of sites and consider the following process for 
revealing it: 

• Initialize S to consist of the sites of T. 

• Order the edges in A arbitrarily {e±, 62, ■ ■ •}• 

• Repeatedly reveal the smallest (as per the above ordering) unexposed 
edge e that is incident to S, setting its value in oj l via an independent 
unit variable U e . If open, add the new endpoint of e to S. 

• Let T be the minimal time after which S has been exhausted (that 
is, the number of edges exposed in the above procedure). 

Upon exposing the sites of S we construct the configurations uj^ ~ i/^ (■ | A^) 
and a/' ~ I A v ) as follows. Let A t denote the set of edges exposed up 
to (including) time t, and suppose that at time t + 1 we are about to expose 
the edge e via the unit variable U e . By definition, 

u\e) = t{U e < v\ ( w (e) = 1 1 w(A t ) = w^At))} , 

and we determine w T (e) for r £ {£, 77} analogously using the same U e : 

oj T (e) = l{U e < i/£(w(e) = 1 1 A T , w(A t ) = w r (A t )) } . 

Beyond time T, we reveal the remaining edges according to their aforemen- 
tioned ordering while following the same recipe. Clearly this construction 
has the correct marginals z^(- | A^),u^{- \ A v ),u^ and the following claim 
establishes that ui l dominates lj^,lj v . 

Claim 3.6. For r G {£, 77} and any integer t the above-defined coupling 
satisfies u) T (A t ) < w 1 (A t ). 

Proof. Let r £ {£,7/}, assume by induction that the statement of the claim 
holds for i and let e be the edge exposed at time t + 1. As this means 
that u; r (Af) < w 1 (At), the Domain Markov property and monotonicity of 
boundary conditions give 

v\ (w(e) = I w(A t ) = w^A*)) < *i (w(e) = | w(A t ) = w T (A t )) . 

Since A 1 " is a decreasing event, the FKG-inequality for the FK-Ising model 
together with the Domain Markov property (imposing a boundary condition 
that arises both from r and from uj T (A t )) ensures that 

v\ (u T (e) = 0,^| W (A t ) = w T (A t )) >^ {u T (e) = | W (A t ) = u T (A t )) 

■^(^|o;(A t )= W T (A t )) , 
and combining the last two equations we deduce that 
4 ( W (e) = I u(A t ) = u \A t )) < v\ ( W (e) = | ^ , w (A t ) = W T (A t )) ■ 
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The statement now follows from the definition of the coupling. ■ 

Observe that if 3 is not incident to the top boundary of A then necessarily 
the edge-boundary between 3 and its complement 3 C is entirely closed in 
cj 1 , and in particular there is a horizontal crossing path in the dual. By the 



monotonicity argued in Claim 3.6 this carries to u T for r G {£, rf\ as well. 



The following RSW-type estimate of 13 for the critical FK- model under 
arbitrary boundary condition will now imply that with positive probability 
3 is confined to the box [0, r + 1] x [0, pr — 1] (as illustrated in Fig. [3]): 

Theorem 3.7 (|13[ Theorem 1]). Let < a.\ < ct2 ■ There exist two 
constants < c~ < c + < 1 ( depending only on a\ and ct2 ) such that for any 
rectangle R with side lengths n and m G [ain, o^ra]] (i-e. with aspect ratio 
bounded away from and oo ) one has 

C < vi(C v {R)) < c+ 

for any boundary conditions £, where C V (R) is the event that there is a 
vertical open path connecting the top and bottom boundaries of R. 

Indeed, applying the above theorem on the box R = [0, r+lj x [0, ar\ with 
wired boundary conditions implies that z^(C„(i?)) < c + for some c + < 1 that 



depends only on the constant a > specified in Theorem 3.2 Introducing 
wired boundary at level ar and applying the Domain Markov property we 
can iterate this argument to obtain that for R = [0, r + 1] x [0, prj we have 
u]i(C v (R)) < (c + )^/ a J with the same c + < 1. Altogether: 

Corollary 3.8. Let r'/r > a > for some a > fixed and let a < p < r'/r. 
There exists some constant 5 = 5(a) such that v%(C v (R)) < exp(— dp) for 
R = [0, r + 1] x [l,pr] and any boundary condition where C V (R) is the 
event that there is a vertical open path connecting the top and bottom of R. 

Since the monotonicity of boundary condition and the Domain Markov 
property imply that vj^ stochastically dominates u^{oj(R) G •), we get 

P (3 c [0, r + 1] x [0, pr - 1]) > 1 - exp(-^) . (3.3) 

Denote by T the edges that are not incident to any of the sites in 3. 
We argue that conditioned on the edges revealed up to time T, namely 
w^(Ay) and oj v (At), the measures z^(-|^) and u^(-\A v ) on the remaining 
unexposed edges T are identical. In other words, 

4(w(T)€-|^,w(A T )=£^(A T )) 

= v\ (w(T) G • | A" , co(A T ) = c^(At)) • (3.4) 

To see this, first consider (w(T) G • | uj(Ax) = oj^(Ax))- Notice that this 
measure over the unexposed edges T is the FK-model on the graph induced 
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on the sites E c , where the edges in the interface with H are all closed and 
the boundary conditions inherited from £ are precisely the pair-wise wirings 
of the sites of V(£) n H c and of those in n H c . Crucially, by our 

construction (that conditioned on A^) there is no open path between V(£) 
and .M (£) in u?(At), hence the sites V(£) D H c are not wired to the same 
cluster of M(£) D H c at time T. 

An analogous statement holds for w^Ay), and since £, rj differ only on T 
we infer that v\ (w(T) E • | u(A T ) = uj t (A t )) has the same boundary con- 
ditions for both r S {£, 77} . The Domain Markov property now gives 

4(^(T) G • I w(A T ) = «f(A r )) = ^(w(T) e • | w(A T ) = w"(A T )) • 

Furthermore, since there was no open path V ~-> .M neither in oj^(At) nor 
in £j ?? (Ar) and the interface between H,H C is closed, the final configuration 
w T for r £ {£, ry} can only have such a path if it is contained in T, which 



together with the last equality implies (3.4) 



Observe that, as a result of (3.4), beyond time T our procedure automati- 
cally couples the configurations io^,uj v via the identity coupling. Altogether, 
we have shown that we can couple u*, uj v on T, in particular obtaining that 
every component C C S c is identical between these two configurations, and 
any component C intersecting C is completely contained in C (since the 
edges between E and H c are all closed). 

We may now assign spins according to Corollary |3.5| Since there is no 



path between V, A4 in either of the configurations, we can assign plus-spins 
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to V and minus-spins to Ai when generating the Ising configurations a^,a v . 
Coupling the i.i.d. spin values for the remaining clusters of S c to be the same 
in both configurations gives <t^(H c ) = a v CE c ). Corollary 3.5 further implies 
that a T ~ pi\ for both r G {£, ??}. 



The proof is concluded by inequality (3.3) which ensures that A T C H c 
with probability at least 1 — exp(— 5p). ■ 

Remark. The exponent C > in Theorem [I] can be readily made explicit in 
terms of the crossing probabilities in the fully-wired critical FK-model. For 
instance, in the setting of the square lattice of side- length n the proof gives 

C = 21og 3/2 [2/(1 , 

where P + {\) is the limiting vertical crossing probability in the FK-model on a 
fully- wired rectangle with conformal modulus ^. Given the limiting function 
p + for all conformal moduli one could improve the resulting exponent by 
optimizing the overlap between the blocks in the recursive analysis. See 



28 , 29 for numerical estimates of the function %r 



4. A POLYNOMIAL LOWER BOUND AND OTHER BOUNDARY CONDITIONS 

In this section we establish a lower bound on the inverse-gap of the 
Glauber dynamics for the critical Ising model on a square lattice and extend 
the upper bound of Theorem [T] to periodic/free boundary conditions. 

4.1. A polynomial lower bound on the inverse-gap. To complement 
our result that the inverse-gap is bounded from above by a polynomial in 
the side-length n we show the following: 

Theorem 4.1. Let gap^ be the spectral-gap of the Glauber dynamics for the 
critical Ising model on a square lattice A of side-length n with an arbitrary 
boundary condition £. Then (gap^) > en 7 / 4 for some absolute c > 0. 
Furthermore, this also holds for rectangles with shorter side-length n. 

Proof. In the static Ising model, as the temperature decreases from high to 
critical, Onsager 144] showed that the correlation of two spins transitions 
from having an exponential decay to a polynomial one. As mentioned in 
the introduction, the dynamics also exhibits such a critical slowdown. At 
criticality, Holley |21| showed that in the infinite volume lattice the spin- 
spin autocorrelation Kao(0)at(0) decays like t -1 / 4 rather than exponentially. 
This can be translated to a polynomial lower bound on gap -1 that is slightly 
worse than linear. Here we will use the spin-spin correlations to obtain the 
n 7 / 4 lower bound. 

Consider the Ising model on A = [l,n]] 2 with some boundary condition 
£ and let A* = [n/4, 3n/4j 2 . Our test-function / will be the magnetization 
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over A*, i.e. f(a) = YlxeA* a ( x )- Recalling the definition of the Dirichlet 
form £ (/) in (1.4), we clearly have £ (/) = 0(n 2 ) since / is 1-Lipschitz. It 



thus remains to estimate Var f ( f ) . 

Analogous to our approach in the previous section, we construct an sample 
of the Ising model at equilibrium by going through the critical FK-model 
via the Edwards-Sokal coupling. We first reveal the set of open bonds in the 
FK-model on A \ A* (that is, every bond with at least one endpoint outside 
of A*) and let F denote the <r-algebra generated by these variables. This 
is followed by revealing the remaining FK-configuration and constructing 
the Ising configuration via the Edwards-Sokal coupling. When revealing the 
FK-configuration we condition on the event A^, i.e. that there is no open 
path connecting plus and minus sites of £. 

By the total-variance formula we have that 

Var(/) = EVar(/ | F) + VarE[/ | F] > EVar(/ | F)1 {A ^ 8A} , 

where {A* <9A} denotes the event that there is no open path connecting 
A* to dA (i.e., there is a circuit of open edges surrounding A* in the dual 
FK-configuration). As observed by |13] , this event has probability at least 
c\ > irrespective of the boundary condition £ without conditioning on A^. 
By the FKG inequality, this conditioning only increases the probability of 
the event {A* ^ dA}. 

Observe that conditioned on the above mentioned event, any site x E A* 
is disconnected from dA and so in the final Ising configuration its expected 
magnetization is Ecr(x) = 0. 



Onsager's decay of correlation result was reobtained by 13 using their 
RSW-type inequalities for the FK-model. In particular, they established 
that if A* is a box of side- length m and x, y G A have distance at least em 
from dA for some fixed e > 0, then for any boundary condition 77 

v\» (x +~* y) > c 2 m- l/A (4.1) 

for some C2 = 02(e) > 0. In our setting this implies that 

Var(/ I F, {A* ^ OA}) = £ E[a(x)a(y) \ F, {A* ^ dA}] 

> ^ ¥(x ^ y\ F, {A* ^ dA}) > c 3 n 15/4 . 

Combining this with the lower bound c\ > on the probability of {A* dA} 
we deduce that Var(/) > C4n 15 / 4 for some absolute C4 > 0. Plugging / as a 
test- function for the spectral-gap in ( 1.3 ) we now deduce that gap^ < cn -7 / 4 
for some absolute c > 0, as required. ■ 
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4.2. Free or periodic boundary conditions. The next theorem estab- 
lishes polynomial mixing for free/periodic boundary conditions. In fact, the 
same method infers this result for any mixed boundary conditions, e.g. pe- 
riodic on one side and a mixture of free and fixed conditions on the other. 
In what follows we restrict our attention to boxes of bounded aspect-ratio 
although our arguments from the previous section can be applied to achieve 



analogues of Theorem 3.1 



Theorem 4.2. Let gap^ be the spectral-gap of the Glauber dynamics for 
the critical Ising model on a square lattice A of side-length n under free 
or periodic boundary conditions. Then (gap^) 1 < n c for some absolute 
constant C > 0. 

Proof. The proof of Theorem [T] holds with slight modifications in the current 
setting and in what follows we describe the required adjustments. 

Consider first the case of free boundary conditions. Through the course 
of the recursive analysis, some of the boundaries of the blocks Ai, A2 may 
be free while others have a fixed (arbitrary) boundary condition (e.g., in the 
first step there are 3 sides with free boundary). Hence, the entire proof holds 



given the following variant of the spatial mixing statement of Theorem 3.2 



Theorem 4.3. The statement of Theorem 3.2 holds when the given bound 



ary conditions £, r\ are possibly free on one or more of the sides of dA. 
The proof of the above theorem proceeds almost exactly the same as that 



of Theorem 3.2, where the only essential difference is that a side with free 
boundary now belongs to A rather than to A. 

Notice that the exponent C > obtained in the above proof for free 
boundary is identical to the one from the proof of Theorem [T] 

It remains to treat periodic boundary conditions. In the Edwards-Sokal 
coupling of the Ising and FK-Ising models, a periodic boundary condition 
corresponds precisely to wiring the two identified boundary sites. Thus, we 
may immediately infer another analogue of the spatial-mixing result special- 
ized to two possibly disagreeing sides under the boundary conditions, with 
the agreeing sides possibly having a periodic boundary condition: 

Theorem 4.4. Let A = [1, r] x [l,r'J for some integers r,r' satisfying 
a\ < r'/r < 02 with 01,02 > fixed and let A M = [l,r] x {ipr' , (1 — if))r'\ 
for some fixed < ifi < |. Let £,77 be two boundary conditions on A that 
differ only on the top and bottom boundaries [l,r] x ({0} U {r' + 1}) and 
possibly have periodic boundary on the remaining sides. Then 

/4HA M ) g •) - n\ (a(A M ) G •) <l-<5, 

TV 

where 5 > is a constant that depends only on a±,a2,ip- 
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To prove the above result, apply the following variant of the argument of 



the proof of Theorem 3.2 on the corresponding FK-Ising model. Instead of 
exposing the set H comprising the sites connected to the bottom boundary 
r by an open path, we simultaneously expose H, H' where S' is its analogue 
with respect to the top boundary. The same reasoning now implies that 
H U E' C A^j with positive probability, allowing the coupling of the two 
FK-measures corresponding to the boundaries £, 77. 

With the above estimate at hand, consider first boxes with one pair of peri- 
odic boundary conditions. Using the block-dynamics with recursive analysis 
as in Theorem [TJ this time we choose each of the two blocks to contain both 
of the periodic boundaries, e.g. when these are the top and bottom take 

Ax = [1, r] x (llr', r'j U [1, ir'J) , A 2 = [1, r] x ([1, fr'fl U [§r', r'fl) . 

After one application of this argument (notice that we did not yet make use 
of the possible periodic boundary conditions in the spatial-mixing result) we 
arrive at rectangles with bounded aspect ratio and fixed boundary condi- 
tions, where Theorem [T] already applies. We deduce that the inverse-gap for 
rectangles with bounded aspect-ratio and one pair of periodic boundaries is 
polynomial in the side-length. 

Finally, to obtain this result for tori (corresponding to periodic boundary 
conditions on all sides), split the torus into two overlapping blocks, each 
with one pair of periodic boundary conditions (in the same manner detailed 
above). At this point, the spatial-mixing result (here applied to periodic 
boundary) asserts that with positive probability the complement of one block 
is coupled within a single step. This in turn implies that the inverse- gap of 
the block-dynamics is uniformly bounded, and the aforementioned bound on 
the single-site dynamics (for boxes with a single pair of periodic boundaries) 
completes the proof. ■ 



4.3. Critical anti-ferromagnetic Ising model. It is well-known that the 
anti-ferromagnetic Ising model on the square lattice is equivalent to the 
ferromagnetic model on the lattice with modified boundary conditions (via 
the transformation that flips the spins at all odd sites, including those in 
the boundary). As our bounds hold for any boundary condition we arrive 
at the following: 

Corollary 4.5. The inverse-gap of the Glauber dynamics for the critical 
anti-ferromagnetic Ising model on the square lattice of side-length n under 
arbitrary boundary condition is polynomial in n. 



MIXING OF CRITICAL ISING ON THE SQUARE LATTICE 



23 



5. Concluding remarks and open problems 

In this work we have established that the inverse-gap of the dynamics 
for the Ising model on the square lattice is polynomial in its side-length, 
with a bound independent of the boundary condition. The proof hinges 
on the recent breakthroughs in the understanding of the critical FK- 
representation for the Ising model. 

Furthermore, we show that on rectangles with different side-lengths 
(whose ratio is not necessarily bounded) , the inverse-gap is bounded by 
a polynomial of its shorter side-length only. 

A natural question that arises from this work is to determine the expo- 
nent of the polynomial growth of the inverse-gap for the critical Ising 
model in 1? . At the present time this remains a formidable challenge. 
Another enticing open problem would be to obtain an upper bound on 
the inverse-gap for the critical Ising model in higher dimensions (Z d for 
d > 3). Here the machinery of SLE is no longer available necessitating 
new ideas for the sought after spatial mixing properties of the model. 
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